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Abstract. We study the properties of the globular cluster 
system of the "old" merger NGC 5018. We detect a rel- 
atively poor globular cluster system, that can be divided 
into at least two populations based on colors: a small pop- 
ulation of blue clusters, aged between several hundred Myr 
to 6 Gyr, the formation of which was probably induced by 
the last interaction in the galaxy group; and one (possibly 
two) population(s) of old(er) clusters that can be associ- 
ated with the original galaxy (ies). 

This composite globular cluster system is as elliptical 
(e = 0.5 ± 0.3) as the best oblate model for the luminous 
component (about E6), but as diffuse as the dark halo 
implied by the same model. In terms of a power law slope, 
we find a fall-off of —1.3 ± 0.4 for the globular cluster 
density profile, while the intensity of the stellar light falls 
off with -2.2 ±0.2. 

We sample the globular cluster luminosity function 
down to V = 24.4 mag, and do not reach the turnover. 
Fits to our data with Gaussians and £5 functions can only 
constrain the turn-over to Vto = 26.3 ± 0.6 mag. The 
total number of globular clusters derived from the lumi- 
nosity function and the surface density profile of the glob- 
ular clusters lies then between 800 ± 270 and 1700 ± 750. 
The amount of young globular clusters is estimated to be 
less than 10% of the total, evidence against a significant 
increase of globular clusters by the merger event in this 
case. 

The specific frequency, for which uncertainties in dis- 
tance and total number cancel, is S = 1.0±0.7, lower than 
expected for a normal elliptical galaxy. However, if the re- 
ported young luminous stellar population would brighten 
NGC 5018 by one magnitude compared with normal ellip- 
ticals, the corrected S value would be about 3, comparable 
with values found for ellipticals in small groups. 
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* Based on data collected at the European Southern Obser- 
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Finally, the study confirms that NGC 5018 contains 
a large amount of dust. The globular clusters might even 
reveal a dust concentration. 

Key words: globular cluster systems - globular clusters 
- elliptical galaxies - galaxies: individual: NGC 5018 



1. Introduction 

With the study of globular clusters in interacting and 
merging galaxies, globular cluster systems gain in impor- 
tance for the understanding of the evolutionary history of 
elliptical galaxies. Properties like the luminosity function, 
specific frequency, morphology, and color distribution of a 
globular cluster system contain information on the globu- 
lar cluster formation history, as well as on the distance and 
history of the parent galaxy. Several hypotheses about the 
formation of globular clusters and globular cluster systems 
are discussed in the literature (e.g. Harris 1991, 1994 and 
Richtler 1995 for reviews, Ashman & Zepf 1992 for forma- 
tion in mergers). Besides primordial formation of globular 
clusters there are also indications that the formation can 
be triggered by interactions between galaxies like merging 
and/or accretion of dwarf galaxies. The study of ellipti- 
cals showing indications of merging events or interactions 
provides the possibility to get a direct insight into the 
building up of globular cluster systems. 

Examples of young globular clusters in recent merger 
galaxies are becoming more numerous (e.g. NGC 3597: 
Lutz 1991; NGC 1275: Holtzman et al. 1992; NGC 7252: 
Whitmore et al. 1993, Schweizer & Seitzer 1993; He2- 
10: Conti & Vacca 1994; NGC 4038/4039, Whitmore & 
Schweizer 1995). NGC 5018, that underwent an interac- 
tion about 1 Gyr ago, might be a similar system, probably 
even older than the previous examples. 

NGC 5018 is the brightest member of a small group 
of 5 galaxies (Gourgoulhon et al. 1992). It is classified as 
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Fig. 1. Image taken from the Digital Sky Survey, showing the 
surroundings of NGC 5018 

a "normal" giant elliptical (E4 in the RSA survey; E3 in 
the RC3 catalog) but shows several indications for a past 
merger event. Some properties of NGC 5018 are summa- 
rized in Table 1. 



Table 1. Properties of NGC 5018, taken from the RC3 (de 
Vaucouleurs et al. 1991), except for (V — I) taken from Goud- 
frooij et al. (1994b) 



Name 


NGC 5018 


RA(2000) 


13 13 01 


DEC(2000) 


-19 31 12 


1, b 


309.9, 43.1 


type 


E3 


V 


10.77 ±0.13 mag 


Vopt 


2794 ± 15 km s" 1 


(B-V) 


0.92 mag 


(V - 7)i nner 


1.05 mag 



Both Carollo & Danziger (1994) and Goudfrooij et 
al. (1994b) investigated in detail the photometric profile 
of NGC 5018 out to ~ 100" and found similar results. 
The color of the galaxy is very blue ((V — I) = 1.05, 
Goudfrooij et al. 1994b), and the color gradient is steeper 
than normally observed in ellipticals. This led Carollo & 
Danziger (1994) to conclude that either the galaxy is very 
blue in the outer parts, or that the whole galaxy is rather 
blue and obscured and reddened in the inner part (where 



the boundary between inner and outer parts is at about 
15" radius). Our study excludes this inner part. Both also 
agree that, despite a probably complicated dust distribu- 
tion in the inner region, the isophotes are disky at all radii 
(note that this is not in contradiction with a merger event, 
e.g. Lima-Neto & Combes 1995). 

NGC 5018 is also peculiar in having one of the lowest 
nuclear magnesium absorption-line strengths for its ab- 
solute magnitude (Schweizer et al. 1990). However, if the 
Mg 2 index profile (Carollo & Danziger 1994), that shows 
a strong decrease down to solar values inside 7" from the 
center, is extrapolated to outer radii it reaches again val- 
ues normal for the luminosity of NGC 5018. This suggests 
the presence of a young stellar population in the central re- 
gion, diluting the Mg 2 index. In that case NGC 5018 could 
be considered as a "classical" merger remnant. This also 
is strengthend by the presence of prominent ripples and 
shells (e.g. Schweizer 1987). Schweizer & Seitzer (1992) 
attribute to NGC 5018 one of the highest fine structure 
parameters, making it one of the best merger candidates, 
and derive a heuristic merger age from U BV colors of 5±2 
Gyr. 

Deep HI observations with the VLA by Guhathakurta 
et al. (1989) show a large HI bridge between NGC 5022 
and MCG 03-34-013, the galaxies situated to the E and 
NW of NGC 5018 (see Fig. 1). The authors claim it to be 
the first direct observational evidence for the formation 
of a shell system when an elliptical galaxy merges with 
a cold disk system, in a (at least) three-body encounter. 
They derive from dynamical arguments a relatively recent 
last interaction, about 600 Myr ago, that produced the HI 
plume, and an orbital time for the outermost shell of 300 
Myr. 

Dynamical investigations show that NGC 5018 has a 
relatively high velocity dispersion (220 km s _ , Davies et 
al. 1987) but a flat velocity dispersion gradient. To re- 
produce this fiat gradient a massive, diffuse dark halo is 
required which is half as flattened (about E3) as the best 
fitting oblate dynamical model of the galaxy light (about 
E6, see Carollo et al. 1995, Carollo & Danziger 1994). Also 
the HI velocity data (Kim et al. 1988) support the ex- 
istence of an substantial amount of dark matter in the 
galaxy group around NGC 5018. 

2. Observations and Reduction 

2.1. The observations 

Deep Bessell V (ESO #451) and Gunn i (ESO #461) 
CCD images were taken in direct imaging mode with the 
Danish 1.54m telescope at the European Southern Obser- 
vatory, La Silla, Chile during the night of 23/24.4.1995 un- 
der photometric and good seeing conditions (between 0"8 
and l'.'l, as measured from the FWHM of stellar images). 
The CCD chip in use was a Tektronix chip (ESO #28) 
1024 2 pixels; the scale is 0.38" per pixel, giving a field size 
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of 6.4' x 6.4'. We obtained respectively five 10 and five 12 
min exposures in V and i centered on NGC 5018. Short 
exposures of 1 min were taken in both filters to compare 
our photometry with the published aperture photometry 
of the galaxy (Poulain & Nieto 1994) and check the zero 
point of our calibration. An additional background field, 
located 11' north of NGC 5018 (see Fig. 1), was observed 
to get an estimation of the contribution of background 
objects in the NGC 5018 field. Bias- and sky-flatficld ex- 
posures were taken at the beginning and end of the night. 
Landolt (1992) standard fields were observed throughout 
the night for the photometric calibration. 

2.2. The reduction 

All the reductions were done in IRAF. Bias subtraction 
and sky flat-fields led to satisfactory (flatter than 1%) im- 
ages. The different long exposures were combined with a 
sigma clipping algorithm to clean the master exposures of 
cosmic rays. 

Object search, photometry, and the determination of 
the completeness factors were done with the IRAF ver- 
sion of the DAOPHOT II package. Before applying the 
final point-source search and PSF fitting, we prepared the 
master exposures in each filter by the following proce- 
dure. First we subtracted from the image the brightest 
stars and extended objects with the IRAF imedit routine. 
Then we computed an isophotal model (using the IRAF 
STSDAS isophote package) that we subtracted from our 
cleaned image. The shells, dust lanes and ripples were then 
visible as residuals. To remove also these structures we 
first subtracted all point-source objects by the standard 
DAOPHOT II PSF fitting. Then we applied a 11 x 11 
pixel median filter to the image. The sum of the isopho- 
tal model and the median image was subtracted from the 
original master exposures, leading to a flat background 
image ready for the final point-source search. We checked 
that no systematic offset was introduced in the photome- 
try by this procedure. Due to model fitting residuals and 
the very complex structure in the center of the galaxy we 
did not consider the innermost 30" radius. 

The calibration was done via typically 40 standard 
stars from the Landolt (1992) list by which we transformed 
our colors to Johnson V and Cousins /. The calibration 
agreed well with the aperture photometry values published 
by Poulain & Nieto (1994). The mean residuals for V and 
/ are 0.014 and 0.015 respectively. Typical errors in color 
at the faintest magnitudes are about 0.19 mag, about 0.14 
mag at V ~ 24.0, and drop below 0.10 mag for V < 23.5. 

The determination of the completeness was done by 
standard artificial star experiments. The completeness 
starts to affect the results at magnitudes of 24.0 mag in V 
and 23.0 mag in I for the NGC 5018 field and 23.8 mag in 
V and 22.6 mag in I for the background field. The limiting 
magnitudes where the completeness factor has dropped 
down to 50 % (we used this limit for our investigations) 



are 24.9 mag in V, 23.7 mag in I for the galaxy field, 24.4 
mag in V and 23.3 mag in I for the background field. 

3. Colors 

In this section we derive the mean colors of the globular 
clusters as well as features in the color distribution. Back- 
ground galaxies can be identified by their second moments 
of intensities down to V = 23.0. In the following we will 
use the term "objects" for objects obeying this coarse se- 
lection, i.e. evident galaxies are rejected. 

The color distribution for all the objects in the field 
detected in V and I is shown in figure 2, where the dashed 
lines mark our 80% completeness limits in the V and / 
filters. 




Field (V-I) Background (V-I) 

Fig. 2. Color magnitude diagram for all objects found in V 
and I in the field centered on NGC 5018 (left panel) and in 
our background field (right panel). Note that the galaxy field 
is deeper by 0.4 mag in V and I than the background. The 
dashed lines mark the 80% completeness limits in V and I. 



3.1. Existence of a dust lane? 

We plotted the spatial distribution of objects of different 
colors. The distribution of blue, (V — I) < 0.8, objects 
as well as red, 1.5 < (V — I) < 2.0, objects is shown in 
figure 3, together with a contour plot of NGC 5018. The 
first group is chosen to be bluer than most background 
objects, and the second to be redder than most globular 
clusters. Two features appear in these plots: blue objects 
are totally missing in a stripe extending from the South- 
East corner to the North of the galaxy (0 objects out of 
57, in a lane representing 18% of the surface), while red 
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Fig. 3. Spatial distribution for objects in the field with 
(V — I) < 0.8 (right upper panel), and 1.5 < (V - I) < 2.0 
(right lower panel), as well as a contour plot of NGC 5018 (left 
upper panel). The stripe extending from the South- East cor- 
ner to the North of the galaxy, showing no blue objects and 
an overabundance of red objects, is indicated in the lower left 
panel. The x and y axis are labeled in pixels of 0.38", i.e. 1000 
pixels are equivalent to 6.3'. 



objects are overabundant in this stripe (34 objects out of 
125, or 27% of the objects on 18% of the surface). 

We checked for any fiatfield or reduction effect, as well 
as any similar feature in the background field, but could 
not find any. We exclude the possibility that the inhomo- 
geneous distribution of objects of different colors could be 
due to random fluctuations in the spatial distribution. The 
probability for a configuration as we see it is very small, 
since the "blue" and "red" distributions shown above rep- 
resent already 3 and 2a fluctuations on their own, and are 
correlated to each other as well as with the HI distribution 
as shown below. We are thus confident that this effect is 
real, and see two possible explanations: 

• different populations of globular clusters have an inho- 
mogeneous distribution in NGC 5018. Or, 

• there is a strong reddening along this stripe, e.g. a dust 
lane. 

The first hypothesis is supported by various facts. As 
will be shown below, an unrealistic assumption for the 
gas-to-dust ratio would have to be made to reconcile the 
HI observations with the existence of a dust lane. More- 
over, the possible extinction could not be clearly seen in 
the galaxy light on our color image (note however, that 
this is not trivial, considering the amount of structure sur- 
rounding NGC 5018 and the fact that it occurs where the 
surface brightness drops). Also, the inhomogeneous spa- 



tial distribution of the objects could be a remnant of the 
merger event: Guhathakurta ct al. (1989) noticed that the 
HI bridge connecting NGC 5022 and MCG 03-34-013 
seems to bifurcate at approximately the position and ra- 
dial velocity of NGC 5018, and while the northern part 
continues unbroken, the southern portion disappears near 
the elliptical. This could have lead to an inhomogeneous 
distribution of the globular clusters, if for example they 
formed from the infalling gas in the southern region, where 
indeed the blue objects are seen. However this would not 
explain the over-density of red objects in the northern 
part. 

The second hypothesis, the presence of a dust con- 
centration, is supported by the following points. The 
large HI bridge between NGC 5022 and MCG 03-34-013 
(Guhathakurta et al. 1989), crosses our frame exactly 
along the stripe suspected to be dust (see their Figure 1). 
However the column density in the bridge is around sev- 
eral 10 19 atoms cm~ 2 (Guhathakurta priv.com.) which, 
using the hydrogen column density to color excess ra- 
tio in our Galaxy of Ahi — 5 x 1Q 21 Eb~v atoms cm~ 2 
(Bohlin, Savage and Drake 1978), would give an redden- 
ing of Eb-v = 0.005, far too small to produce the effect 
of several tenths of magnitudes in Ey-i observed, unless 
the gas-to-dust ratio is higher by about two orders of mag- 
nitude (i.e. Mg as /Mdust =1 to 5) than in our Milky Way. 

Further, emission at A60 ^m and A100 /xm at an in- 
tensity of about 1 Jy was detected with the IRAS satel- 
lite. The derived temperature of the infrared emission was 
below 30 K, leading to the assumption that this ellipti- 
cal galaxy possibly contains cold interstellar gas (Knapp 
ct al. 1989). Unfortunately, even the HIRAS maps (pro- 
cessed IRAS images for a higher resolution of 1 arcmin) 
at A60 yum and A100 ^m do still not have the required 
resolution to locate the dust in NGC 5018. 

Finally, the isophotal shape and the distribution of ion- 
ized gas and dust have been extensively studied by Goud- 
frooij et al. (1994a, 1994b) and Goudfrooij & de Jong 
(1995). They noticed that the overall visible optical ex- 
tinction is about 1 magnitude smaller than expected from 
the high IRAS flux, and might be caused by a diffuse dis- 
tribution of dust, probably accreted from other galaxies. 

We have at this stage no reason to prefer one scenario 
to the other and will have to wait for ISO observations in 
order to solve the ambiguity. Nevertheless, we will show 
that most conclusions are not severely affected by the as- 
sumption of one or other scenario, but that in details our 
results might be affected by unknown extinction being 
clumped or diffuse. 

3.2. Color distributions 

We present in the following the color distributions for both 
hypotheses. For the first, we left all colors unchanged. For 
the second, we made a coarse correction for reddening 
within the stripe shown in Fig. 3. 
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3.2.1. The sample 

To determine the color of globular clusters, we selected 
a sample of objects in which the globular clusters domi- 
nate the background and compared this sample with an 
equally selected one on the background field. This sample 
consists of all objects in a elliptical ring from 30" to 170" 
major-axis, the lower limit protecting against variation in 
completeness towards the center, and the upper limit be- 
ing chosen from the density profile: out to a semi-major 
axis of 170", globular clusters dominate the background 
(see Sect. 5.1). 

Further, only objects fainter than V = 22.5 mag are 
considered. This is motivated by the globular cluster lumi- 
nosity function (see Sect. 4.1), which shows that globular 
clusters are only to be expected at fainter magnitudes. 

3.2.2. Correction for eventual local extinction 

For the color correction in the stripe we assumed that 
the colors outside the stripe are unaffected. We corrected 
the colors within the stripe by a constant value to match 
their mean with the mean color outside the stripe. In the 
stripe, the objects following the selection described above 
have a mean color (V — I) = 1.68 mag ; in the rest of the 
ring, the mean color is (V — I) = 1.38 mag. Wc therefore 
adopted a "reddening correction" of A(V — /) = 0.3 mag, 
AV = 0.55 mag , and AI — 0.25 mag. In the following, 
colors and magnitudes in the stripe will be corrected by 
these values when we mention a correction for the "dust 
lane". 

3.2.3. The color distribution 

We plotted the color distribution for objects from V = 
22.5 to 24.5 mag in Figs. 4 and 5. Figure 2 shows that 
completeness corrections start to be necessary at magni- 
tudes fainter than V — 23.5, and that the correction is a 
function of the colors of the objects. 

We divided our sample into four categories: (V — I) = 
0.5 to 1.0, 1.0 to 1.5, 1.5 to 2.0, and 2.0 to 2.5. We then 
applied the correction as a function of magnitude and of 
color to the bins in a given category, computed as the 
product of the mean V completeness and the correspond- 
ing mean / completeness^. See Tab. 2 for the values used 
for the field and the background. The color distribution of 
the background was smoothed by the median over respec- 
tively 3 bins to avoid any large fluctuations due to galaxy 
clusters that we could identify on our background frames, 
which fall in single color bins. It was then scaled down to 

1 We examined in artificial star experiments the effects of 
combining completeness factors for V and I. In our case the 
completeness drops due to photon noise rather than to crowd- 
ing, so that completeness in V and / are independent. We 
therefore used the product of the completeness factors in V 
and in I, rather than just the lowest completeness at a given 
magnitude. 



Table 2. Completeness values in percent as a function of mag- 
nitude and color for the field (left value) and the background 
(right value). 



Mag. / (V -I) 


0.5—1.0 


1.0—1.5 


1.5—2.0 


2.0—2.5 


23.0-23.5 


95 95 


100 100 


100 100 


100 100 


23.5-24.0 


87 77 


95 95 


100 100 


100 100 


24.0-24.5 


68 22 


87 54 


88 66 


91 68 



the ring area considered in the galaxy field and subtracted 
bin by bin from the color distribution of objects in the ring 
around NGC 5018. 




0.5 1 1.5 2 2.5 



(v-i) 

Fig. 4. Color distribution of globular cluster candidates be- 
tween 30"and 170", with magnitudes 22.5 < V < 24.5. The his- 
togram shows the final corrected counts, the histogram of the 
background used for correction is over-plotted (dashed line). 
No correction for the "dust" lane was applied. 

The resulting color distribution is much broader than 
expected from the photometric errors alone. If all clusters 
were concentrated at a single color, we would obtain a 
gaussian distribution with a FWHM of about 0.4 mag. 
The distribution uncorrected for reddening shows different 
groups of objects: a small group of objects around (V — 
I) = 0.7 mag, a broad population around (V — I) — 1.2 
mag, and a narrower one around (V — I) = 1.65. 

The color distribution corrected for the reddening in 
the "dust lane" shows a different mixing of the popu- 
lations. The population at (V — I) = 0.7 becomes as 
large as the one at (V — I) = 1.2, and the population 
at (V — I) = 1.65, shifting to a somewhat bluer color 
((V — I) = 1.50), is only seen as a tail of the population 
at (V-I) = 1.2. 
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Fig. 5. Color distribution of excess objects between 30"and 
170", with magnitudes 22.5 < V < 24.5. The histogram shows 
the final corrected counts, the histogram of the background 
used for correction is over-plotted (dashed line). The colors in 
the stripe (see figure 3) were corrected by A(V — I) = 0.3 mag. 

Both distributions show some red objects clustered 
around (V — 7) = 2.2 or (V — I) = 1.9 respectively. In 
the case of no reddening, they would be a galaxy clus- 
ter at low redshift, judging from their color. In case of a 
reddening stripe the whole group shifts towards the blue: 
all these objects appear to lie in the stripe. This shows 
that our reddening correction of 0.3 mag might be very 
crude, and that, in fact, some objects could be even more 
reddened. 

3. 3. Interpretation of the color distribution 
3.3.1. The colors 

In the following we include the Galactic extinction of 
E(B -V)= 0.053 mag or E(V - I) = 0.072 mag towards 
NGC 5018 (Burstein & Heilcs 1984). 

If we compare the colors of the globular clusters in 
NGC 5018 with those of the Galactic globular clusters (see 
Fig. 6, data taken from Harris 1996), we notice that the 
distribution of the two red populations is similar to that of 
the globular clusters in the Milky Way, if the colors of the 
latter are not corrected for extinction. From their colors, 
the red population(s) are thus compatible with being old 
globular clusters as seen in the Milky Way, with metallic- 
ities ranging from [Fe/H] = —2.5 to dex and possibly 
suffering from as much reddening. We cannot discriminate 
between single populations with a wide range of metallic- 
ity combined with some reddening due to the dust in NGC 



i 1 1 1 1 i 1 1 1 1 i 1 1 1 1 i 1 1 1 r 




(v-i) 



Fig. 6. Color distribution of the Galactic globular clusters, un- 
corrected for reddening. The finer hashed part is the contribu- 
tion of the globular clusters for which the reddening E(B — V) 
is less than 0.2. Data from Harris (1996). 



5018, or eventually two distinct populations of old globu- 
lar clusters. 

However, the blue objects are bluer than even the 
most metal-poor globular clusters in the Milky Way. 
Therefore they must be younger, unless they are more 
metal-poor than any observed globular cluster to date 
(e.g. Geisler et al. 1995). According to the models from 
Fritze v.Alvensleben & Burkert (1995), agreeing with 
Bruzual & Chariot's (1993) model for solar metallicity, 
the globular cluster population centered at (V — 7)o = 
0.63±0.10 (eventually redder by 0.1 in the corrected case), 
could have an age ranging from 200 Myr to 6 Gyr, depend- 
ing on the metallicity of the globular clusters (the younger 
limit being for extremely metal rich (Z — 4 x 10~ 2 ), the 
upper for extremely metal poor [Z = 1 x 10~ 4 ) globu- 
lar clusters). This compares well with the young globular 
clusters found in NGC 7252 (Whitmore et al. 1993), which 
are marginally bluer (by about 0.1 mag) and brighter by 
2-3 magnitudes than our candidates. And thus probably 
younger by a few 100 Myr. These clusters are estimated 
to have ages between and 800 Myr (Schweizer & Seitzer 
1993). 

All this supports the fact, that the young globular clus- 
ters in NGC 5018 could have formed during the last inter- 
action dated between several hundred Myr and a few Gyr 
ago (see references in Sect. 1). 



Michael Hilker & Markus Kissler-Patig: The globular cluster system of an "old" merger: NGC 5018 



7 



3.3.2. The color distribution with magnitude 

In Table 3 we show the counts of globular clusters (back- 
ground subtracted) for two magnitude bins, in order to 
examine the distribution of the red and blue objetcs with 
the magnitude. For the two red intervals the counts in- 



Table 3. Number of globular clusters after background sub- 
traction in different color intervals (see Figs. 4 to 5 for the 
histogramms) for the case where no dust lane is taken into ac- 
count, and for the case where coorections are made for a dust 
lane as explained in Sect. 3.2.2. 



color interval: 


0.6—0.9 0.9—1.4 


1.4—1.9 


mag limit 








uncorrected case 




22.5—23.5 


1.4 13.6 


13.2 


22.5—24.5 


13.2 56.8 


42.3 




corrected case 




22.5—23.5 


4.5 12.7 


7.6 


22.5—24.5 


28.8 52.1 


29.7 



crease by a factor 4 when going a magnitude deeper (in- 
dependently of a correction or not) . Assuming a turnover 
magnitude around m v — 26.3 mag (see Sect. 4.1), we cut 
off respectively 2.8 and 1.8 mag below the turnover mag- 
nitude and sample respectively 3% and 12% of the surface 
of the globular cluster luminosity function (assuming it to 
be Gaussian and to have a dispersion of 1.4 mag). This is 
in good agreement with the increase by roughly a factor 
4 found for the red populations. 

For the blue objects, the statistics are too uncertain to 
be used in a quantitative way. Qualitatively, we detect only 
a few (1-5) objects down to V = 23.5 (including also the 
interval V = 20.5 to 22.5, i.e. M v > -13, adds 3 more ob- 
jects) and then a sudden increase (13-30) when sampling 
a magnitude deeper. Assuming an unimodal luminosity 
function for the blue objects, not much broader than the 
luminosity function of red objects, this would hint at a 
peak about 2 3 magnitudes brighter in V than that of 
the red objects. This fits into the scenario where the blue 
objects would indeed be younger: according to e.g. Fritzc 
v.Alvcnsleben & Burkert's model, globular cluster around 
1 Gyr would be 3 magnitudes brighter in V than 15 Gyr 
old ones. 

To get an estimate of the ratio of blue to red globu- 
lar cluster candidates we assume a "blue GCLF" with a 
Gaussian shape, peaking 2 to 3 magnitudes brighter than 
the GCLF of the red objects, for which we take the turn- 
over values from table 4. This leads then to 20-150 blue 
globular clusters in total, or 1-10% of the red objects (see 
Sect. 4.2). 

In summary, despite the fact that we have small num- 
ber statistics for the globular clusters and large uncertain- 



ties in the dust distribution, the following points can be 
retained: 

• A blue population of globular clusters exists in NGC 
5018, and is probably between 200 Myr to 6 Gyr old. (the 
upper range would be slightly preferred from their mag- 
nitudes, when compared to the young clusters in NGC 
7252). These young globular clusters are likely to have 
formed in the last interactions within the galaxy group. 

• The red globular cluster population(s) is compatible 
with old (> 10 Gyr) globular clusters. Two distinct popu- 
lations might exist, but the data are still compatible with a 
single population differentially reddened and with a large 
spread in metallicity. 

4. The number of globular clusters 

The globular cluster luminosity function (GCLF) plays an 
important role as extragalactic distance indicator and is 
essential for the estimation of the total number of globular 
clusters and specific frequencies. In recent years a lot of 
work has been done towards understanding the represen- 
tation of the GCLFs and the fine tuning with the depen- 
dences on metallicity and age (e.g. Seeker 1992, Ashman 
et al. 1995). 

4-1. The globular cluster luminosity function 

To compute the GCLF, we used the coarse selection of 
the objects described in the Sect. 3. (i.e. obvious galaxies 
were rejected) and restricted ourselves to an elliptical ring 
from 53" to 170" major-axis. We made no selection in 
color, since we wanted to avoid a combined completeness 
correction. Thus we included the blue objects and one has 
to be aware that these objects may influence the following 
results as discussed later. 

To determine the GCLF, we counted the selected ob- 
jects in bins of 0.5 mag width and shifted the bins in steps 
of 0.125 mag to get four luminosity distributions. We then 
corrected the counts for incompleteness. These steps were 
done for the field centered on NGC 5018 as well as for 
the background field. We interpolated the background lu- 
minosity function by a polynomial fit of the 4th order 
to smooth out the effect of individual background galaxy 
clusters biasing given bins. This background luminosity 
function was then subtracted from the field luminosity 
function bin by bin. 

The limiting magnitude of 24.4 mag in V is determined 
by the last bin of the background luminosity function with 
a completeness above 50%. Since we only have 7 to 8 data 
points for each fit to the bright tail of the luminosity func- 
tion, we combined respectively the two luminosity func- 
tions shifted by 0.25 mag to obtain distributions with 0.5 
mag bins in 0.25 mag steps. The points are then no longer 
statistically independent, but stabilize the fit. 

The best representation of GCLFs are Gaussians or 
Student t 5 functions with dispersions of around txcauss = 



8 



Michael Hilker & Markus Kissler-Patig: The globular cluster system of an "old" merger: NGC 5018 



1.4 or crt 5 = 1.1 respectively (e.g. Seeker 1992). The re- 
sulting turnover values (error-weighted means) from fit- 
ting these functions to the double-binned distributions 
vary a lot as a function of the chosen dispersion, while 
being almost unaffected by the correction for the dust 
lane or different binnings. In case of a Gaussian fit the 
values range from Vto — 25.9 for the lowest dispersion of 
CGauss = 1-2 to Vto = 26.9 mag for the highest dispersion 
of 

CGauss — 1-5. For a fit with a t$ function the values 
range from from Vto = 25.7 for the lowest dispersion of 
(Tt s = 1.0 to Vto = 26.4 mag for the highest dispersion 
of tJt 5 = 1.3. The turnover values have been corrected for 
a foreground extinction of Ay = 0.17 mag (Burstein & 
Heiles 1984). 

The uncertainties in our results are the following: 

- The accuracy of the turnover determination by a given 
function lies between 0.3 and 0.6 mag for a sparse sam- 
ple like ours, covering only the bright end of the GCLF. 
This was derived from Monte Carlo simulations, simulat- 
ing samples similar to ours and fitting them in the same 
way. Single deviating bins can largely influence the fit. 

- Figure 7 shows that the fits of a Gaussian and a t$ func- 
tion with different dispersions are indistinguishable from 
each other in the fitted bright tail, whereas the turnovers 
cover a range of about one magnitude. Since we do not 
know the real dispersion for the GCLF of NGC 5018, we 
can only derive turnover values within this range. Never- 
theless, in the case of a broad color distribution one would 
expect a dispersion around OGauss = 1-4, as seen for ex- 
ample in M87 which has a comparable color distribution 
(Whitmore & Schweizer, 1995). 

- A neglected young globular cluster population could af- 
fect the bright tail of the luminosity function. The possible 
young population which is 2-3 mag brighter than the old 
one would produce a more extended, flatter wing on the 
bright side of the GCLF which leads to underestimation 
of the turnover values for the old population. Monte Carlo 
simulations showed that in our case (less than 10% young 
clusters and only fitting the bright tail) the underestima- 
tion would be less than 0.1 mag, but would introduce an 
additional scatter of 0.5 mag. 

Including these uncertainties we can estimate a dis- 
tance for NGC 5018 using the GCLF. We assume an ab- 
solute turnover value corrected for metallicity of M v = 
—7.1 ± 0.2 mag (following Ashman et al. 1995, combined 
with a metallicity slightly less than solar for our old pop- 
ulation, assumed from our colors). 

Our average values lead to a distance modulus in the 
range of (to - Mf v = 32.8 ± 0.4 to (to - M) v = 34.0 ± 0.6, 
bracketing the value given in Tully (1988), derived from 
the heliocentric velocity of v = 2794 ± 15 km s" 1 , and 
assuming H = 75 km s _1 Mpc -1 . 

An accurate distance determination by the GCLF 
awaits deeper data. 
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Fig. 7. The globular cluster luminosity function of NGC 5018 
fitted by Gaussians and t 5 functions with different dispersions 
(see Sect. 4.1). In the bright tail all of the fits are equally good 

4-2. The total number 

The calculation of the total number of globular clusters 
strongly depends on the fraction of the GCLF that we 
have covered with our data. The fitted luminosity function 
covers all clusters in the ring of 53" to 170" of major-axis. 
The results of the Gaussian fits of Tab. 4 (with "dust" 
lane stripe) range between 500 ± 200 and 1000 ± 600 ob- 
jects. For the inner and outer parts one has to make as- 
sumptions for the number of globular clusters based on 
the radial distribution of the globular cluster density (see 
section 5.1). Due to disruption processes we do not expect 
a steep increase of the cluster density in the innermost 10 
kpc from the center but a flattening of the density pro- 
file (e.g. McLaughlin et al. 1993). Therefore, a reasonable 
guess for the inner ring (0" - 53") is a density of 29.3 ± 5.0 
objects per square arcmin (the density at 40") leading to 
50 ± 9 objects down to a magnitude of V = 25 mag. For 
the outer part (170" - 200") we assumed a mean density 
of 4.0 ± 1.7 leading to 27 ± 12 objects. We then corrected 
for the full GCLF. Our estimate of the total number of 
globular clusters in NGC 5018 ranges from 800 ± 270 to 
1700 ± 750. 

4-3. The specific frequency 

The specific frequency is defined as S — N x io°- 4 ( M "+ 15 ), 
where N is the total number of globular clusters in the 
galaxy and My its absolute visual brightness (Harris & 
van den Bergh 1981). Mean values for S are found to 
correlate with the galaxy type, luminosity, and environ- 
ment. While the mean specific frequency of spirals is 
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around 1, ellipticals in small groups have a higher mean 
value around 2.6, and ellipticals in clusters have a mean 
value around 5.4 (Harris 1991). Kumai, Hashi, and Fuji- 
moto (1993) tried to quantify empirically the correlation 
of the specific frequency S with the local galaxy density. 
For the given galaxy density around NGC 5018 of 0.29 
galaxies/Mpc 3 (Tully 1988), their relation would lead to 
a value of S = 3.0 ±0.2, also expected from Harris' (1991) 
compilation. Djorgovski & Santiago (1992), supported by 
Zepf, Geisler & Ashman (1993) showed the relation be- 
tween galaxy luminosity and number of globular cluster 
to be N gc oc L 2 rather than N gc oc L 1 , from which we also 
would expect an S value around 3 for NGC 5018. 

The absolute visual brightness of NGC 5018 depends 
on its distance, but since the derived total number of 
clusters, taken from the GCLF, also increases with dis- 
tance, the calculated S values lie in a narrow range be- 
tween S = 0.9 ± 0.5 and S = 1.1 ± 0.6, significantly below 
the expected value of S = 3.0 for our given N and My 
values. For NGC 5018 to fall in the mean for ellipticals 
in small groups this would mean: cither having 3 times 
as many clusters as observed, or being roughly one mag- 
nitude "too bright". One explanation could be that the 
absolute luminosity of NGC 5018 is increased by a young 
stellar population and therefore S is underestimated. 

We calculated the increase of the luminosity and the 
change in color when mixing a young blue and red old pop- 
ulation based on the models given by Fritze v.Alvensleben 
& Burkert (1995). If one wants to reach a S value of 3 the 
young population would have to have increased the total 

V luminosity of NGC 5018 by about 1.1 mag. A 3 Gyr old 
population with a color of V — I = 0.75 mag which con- 
tributes 33% to the total V luminosity would increase the 
luminosity by that amount. The color would change from 

V — I = 1.25, which is usual for ellipticals (e.g. Goudfrooij 
et al. 1994b), to V - I = 1.05, the color of NGC 5018. 

5. Spatial distributions 

5.1. The radial distributions 

In the following we derive the radial density profile of the 
globular clusters and compare it with the profile of the 
galaxy light. 

We compute the radial profile using all the objects 
found on our frame down to a magnitude of V = 25.0. 
The objects were binned into elliptical rings 100 pixels 
(37.7") wide, the counts were then corrected for finding 
incompleteness (using 100% down to V = 24.0, 94% be- 
tween V = 24.0 and 24.5, and 65% between V = 24.5 
and 25.0), and geometrical completeness when the ring 
expanded over the limits of our frame (only necessary for 
the ring 500 to 600 pixels (189" to 226")). 

Table 4 shows the mean semi-major axis of the ring, 
the corrected counts, the applied geometrical complete- 
ness, the total ring area, and the density. The density error 



was taken as a combination of the square root of the raw 
counts and uncertainties in the correction for incomplete- 
ness. The second part of Table 4 shows the object densities 
on the background field, where no rings were denned, but 
3 stripes perpendicular to the minor axis. 



Table 4. Radial counts of objects around NGC 5018 (upper 
part) and in the background field (lower part) 



radius 


corrected 


geom. 


full area 


density 


in arcsec 


counts 


corr. 


in arcmin 2 


obj /arcmin 2 


56.5 


105.6 


1.00 


2.61 


40.5 ± 3.9 


94.3 


121.4 


1.00 


4.34 


28.0 ±2.5 


132.0 


148.0 


1.00 


6.08 


24.3 ±2.0 


169.7 


200.2 


0.96 


7.82 


25.6 ± 1.8 


207.4 


182.4 


0.79 


9.56 


19.1 ± 1.6 


290.3 


280 


1.0 


14.7 


19.0 ± 2.5 


444.9 


250 


1.0 


13.3 


18.8 ±2.5 


662.8 


254 


1.0 


13.3 


19.1 ±2.5 



On the background field, the counts at V — 24.5 mag 
are already incomplete by more then 50%. We therefore 
corrected these counts by comparing them with Smail's et 
al. (1995) deep counts in V, and extrapolated our results 
down to V = 25.0 mag with their slope. We find for our 
background field a mean density of 19.0 ± 1.2 objects per 
square arcmin down to V = 25.0 mag, in agreement with 
Smail's et al. (1995) results. 

The radial profile is plotted in Fig. 8. We fitted a 
power-law of the kind p <~ r~ x , where p stands for the 
surface density of the globular clusters or the light flux of 
the galaxy, and r for the semi-major axis, to both globular 
cluster and light profiles in order to compare their fall-off. 

For the density profile of the globular clusters after 
subtraction of the background density, we obtain x = 1.3± 
0.4 (the large error being due to the weight attributed to 
the fourth point). The fit to the light profile (between 10" 
and 160" semi-major axis) returns a slope of x = 2.4±0.1, 
and is therefore steeper than the globular cluster profile. 
This last slope depends slightly on how far one extends 
the fits, since the galaxy light gets steeper beyond a semi- 
major axis of about 100". Stopping our fit there returns 
x = 2.2 ±0.1, essentially the same value that would return 
a fit to the Goudfrooij et al. (1994a) data in this range, 
for which we find x = 2.15 ± 0.05. 

The galaxy light does not have an unexpectedly steep 
profile, it is rather the globular cluster system which ex- 
tends further out and is flatter than the stellar compo- 
nent. A similar behavior was already observed in galaxies 
with rich globular cluster systems, while in sparser sys- 
tems the globular clusters seem to follow the galaxy light 
(e.g. Kissler-Patig et al. 1996 and references therein). 
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Fig. 8. Radial density profile of the globular clusters. The up- 
per panel shows the raw object densities in counts/arcmin 2 
down to V = 25.0 mag (circles are from the field, triangles 
from the background) ; the dashed line marks our mean object 
density on the background field. The lower panel shows the 
background corrected densities, i.e. the density of the globu- 
lar clusters; here we over-plotted with an arbitrary offset the 
galaxy light profile (squares). 

5.2. Ellipticity and position angle of the globular cluster 
system 

Figure 9 shows the angular distribution of all objects down 
to a magnitude of V = 25 around NGC 5018 in a ring from 
38" to 188", uncorrected for finding incompleteness; that 
we assumed to be homogeneous over the position angle. 
Objects closer than 38" were not considered because of 
the too high finding incompleteness within this radius. 
To avoid any geometrical incompleteness, wc limited the 
counts to 188". The thus defined ring area covers 29.6 
arcmin 2 , or 73% of the total frame surface. 

We divided the ring into 12 segments, 30 degrees wide 
each, and counted the objects in the segments. Each ring 
has a surface of 2.48 arcmin 2 . We estimated the back- 
ground contamination, uncorrected for completeness, from 
the area outside 190", where the density profile reaches al- 
ready the background (see Sect. 5.1.) and obtained 12.0 
objects per square arcmin. We thus expect 30 objects in 
each segment (marked as a dashed line in Fig. 9). 

The two strong peaks in Fig. 9 separated by 180 de- 
grees at the position angle of the galaxy hint at a non- 
spherical distribution of the objects. Indeed, the distribu- 
tion can be better fit by a double cosine (expected for an 
elliptical distribution) than by a constant function (ex- 
pected for a purely spherical distribution): note however 



Fig. 9. Angular distribution of the objects around NGC 5018. 
The dashed line shows the expected background level, the solid 
curve is the best fit of a double cosine to the data. 

that a spherical distribution is still compatible with the 
data within 2a. 

The best fit by a double cosine of the form n = A * 
cos(2 *((p — B)) + C, where n are the corrected counts and 
4> the position angle E of N, returns the parameters A = 
4 ±2, B = 114 ±25, and C = 11 ±2. The globular cluster 
system thus has a position angle of 114±25 degrees E of N, 
and its ellipticity as 1 - (C - A)/(C + A) = 0.5 ± 0.3. The 
isophotcs of the galaxy have a position angle of 94 ± 2 out 
to 30" and then become more irregular, shifting between 
95 and 110 degrees E of N (Goudfrooij et al. 1994b). The 
ellipticity of the isophotes is roughly 0.3 ±0.1 (Goudfrooij 
et al. 1994b), while the best oblate model for the galaxy 
light returns 0.6 (Carollo & Danziger 1994). 

The globular cluster distribution is compatible within 
the errors with the galaxy light. It seems slightly shifted 
with respect to the stellar light and would rather support 
Carollo & Danziger's (1994) model where the galaxy is 
intrinsically flat (E4 to E6) but seen at an inclination of 
30 to 55 degrees, if we assume that the globular clusters 
share the ellipticity of the luminous component. 

6. Discussion and Conclusions 

6.1. Verifying the prediction of globular clusters in merg- 
ers 

At least some ellipticals are supposed to have formed by 
mergers or interactions. When the majority of these events 
happened is still unclear, as well as the evolutionary state 
of the progenitors (still gaseous proto-galaxies or evolved 
stellar systems?, see e.g. White (1994) for a review). 
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NGC 5018 is probably one of the elliptical galaxies that 
were formed by mergers or interactions of evolved stellar 
systems. Guhathakurta et al. (1989) suggest an encounter 
of at least three bodies including at least one cold disk 
system. 

The properties of the globular cluster system of NGC 
5018 can be compared with the predictions of the model 
for globular cluster formation in mergers (Ashman & Zcpf 
1992 and Zepf & Ashman 1993). Several predictions of the 
models are confirmed: 

(i) - NGC 5018 most probably contains "young" (sev- 
eral hundred Myr to 6 Gyr) globular clusters that could 
be attributed to the last interaction of NGC 5018. This 
would close the gap between different "old" populations 
in the same galaxy and "newborn" globular clusters, as 
seen directly in merging or interacting galaxies (e.g. NGC 
3597: Lutz 1991; NGC 1275: Holtzman et al. 1992; NGC 
7252: Schweizer & Seitzer 1993; Hc2-10: Conti & Vacca 
1994). 

(ii) In NGC 5018, the problem caused by 
age/metallicity ambiguity in broad-band filters is partly 
solved by the big differences in age between the clusters, 
that make the young population clearly stand out in color. 
It remains however very difficult to decide whether one 
or two "old" populations are present, owing to the large 
uncertainties in the dust distribution within the galaxy. 
Therefore it is not possible to say whether two similar 
systems (e.g. spiral galaxies) merged, in which case differ- 
ent old populations are expected, or that a larger system 
(already elliptical?) interacted with smaller systems. 

Different populations of globular clusters around el- 
liptical galaxies, as identified in color distributions (e.g. 
examples in Zepf et al. 1995, or the recent study by Whit- 
more & Schweizer 1995 in M87), wait for confirmation by 
spectroscopy of whole globular cluster systems with 10m- 
class telescopes. 

(iii) - The model also predicts flatter radial density 
profiles for globular cluster after a merger event. The as- 
sumption is that old and new globular clusters follow the 
distribution of the field stars with which they formed. Dur- 
ing a merging event the gas condenses to the center to form 
the new stars and clusters, so that the resulting stellar 
distribution is more centrally concentrated than the old 
globular clusters. If, after the interaction, the old globular 
clusters still dominate the system, it will appear flatter 
then the galaxy light. If, on the contrary, the new globu- 
lar clusters dominate, the system will only be marginally 
flatter then the light of the galaxy. 

In the case of NGC 5018, the globular cluster system 
extends further than the stellar component. Therefore, if 
the assumption above is right, we would expect only very 
few globular clusters to have formed newly and the old 
globular clusters to still dominate the system, as is in- 
deed the case (see Sect. 4.1), since we estimate the young 
population to be at most 10% of the old one. 



In that respect, however, poor systems that follow the 
profile of the galaxy light like NGC 720 (Kissler-Patig et 
al. 1996) would not fit into a merger scenario. Density pro- 
files of globular clusters could then serve to identify ellip- 
ticals that did not form by late stellar mergers. 

One prediction of the model, however, does not seem 
to be verified here: originally, the model aimed to explain 
the higher specific frequency of ellipticals compared with 
spiral galaxies. Taking into account the different mass to 
V luminosity along the Hubble sequence, the specific fre- 
quencies of ellipticals are higher by a factor two than in 
spirals (Ashman & Zepf 1992). If this excess of globular 
clusters had to be explained by formation during mergers 
and interactions alone, Ashman & Zepf (1992) expected 
the formation of as many new globular clusters as there 
were old ones. 

The specific frequency of NGC 5018 is not much higher 
than in spirals, mainly due to its high luminosity. To rec- 
oncile its S- value with that of normal ellipticals in groups 
(about 3, Harris 1991) the young stellar population must 
have raised the absolute luminosity of NGC 5018 by one 
magnitude. 

The results for NGC 5018 show, however, that only 
few blue globular clusters formed in this last interaction. 
We estimate the young population of globular clusters to 
make up less than 10% of the total population. At least in 
this case, the merger event did not produce a sufficiently 
high number of new clusters to systematically explain the 
higher number of globular clusters found in ellipticals as 
compared to spirals. Either NGC 5018 experienced more 
mergers in the past that were more effective, or another 
scenario is needed to explain the systematic higher number 
of globular clusters around ellipticals. 

Harris (1995) estimated that for an efficiency of 1% 
to turn gas into bound stellar clusters during a merger, 
typical merging events involving 10 10 M of gas would 
produce about 300 clusters with an "average" mass of 
3 x 10 5 M Q . In the case of NGC 5018 roughly 100 new 
clusters formed, which might allow an estimation of the 
amount of gas involved. Note that the efficiency seem to 
be influenced by the amount of gas in both merging galax- 
ies. NGC 5018, probably an elliptical merging with a cold 
disk system (Guhathakurta et al. 1989), resembles in that 
respect NGC 5128, the product of a bulge dominated sys- 
tem with a late-type spiral (Minniti et al. 1996), for which 
Zcpf & Ashman (1993) suggested that the small amount 
of newly formed clusters was due to the swallowing of a 
small satellite. In contrast, systems like NGC 4038/4039 
(the "Antennae"), the result of a Sb and a Sc galaxy merg- 
ing, produced over 700 young globular clusters (Whitmore 
& Schweizer 1995). This suggests that "gas on gas" in- 
fall creates very favorable conditions for the formation of 
globular clusters, while the infall of gas into an early-type 
galaxy is not as effective. 
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6.2. Globular clusters compared to dark matter 

In their kinematical study, Carollo & Danziger (1994) 
found through dynamical modeling that NGC 5018 is 
much better fit by a model with a dark matter halo than 
without. Their best oblate model with dark matter re- 
turns the presence of a dark halo six times as massive, 
twice as diffuse, and half as flattened (about E3) as the 
luminous galaxy. In that respect it is interesting to note, 
that the profile for the globular cluster system was com- 
puted further than 38" (> 2R e ), where the dark matter 
potential will dominate, and is also twice as diffuse as the 
stellar component. However, its ellipticity seems higher 
than that returned by the dynamical model for the dark 
halo, even if both still agree within the errors. 

The association of globular cluster systems with the 
stellar or the dark halo also waits verification from kine- 
matical studies of entire globular cluster systems. 

Acknowledgements. We wish to thank Tom Richtler for his 
constant support, as well as Paul Goudfrooij and Edwin 
Huizinga ("the unbelievers") for helpful discussions. Uta 
Fritze-v.Alvensleben is thanked for providing us with the tab- 
ulated results of her globular cluster color evolution computa- 
tions. Thanks also to Do Kester for providing the HIRAS im- 
ages of NGC 5018. The referee S. Zepf is acknowledged for use- 
ful comments. This research was partially supported through 
the Deutsche Forschungsgemeinschaft under grant Ri 418/5- 
1. This research made use of the NASA/IPAC extragalactic 
database (NED) which is operated by the Jet Propulsion Lab- 
oratory, Caltech, under contract with the National Aeronautics 
and Space Administration. 

References 

Ashman K., Zepf S.E., 1992, ApJ 384, 50 

Ashman K.M., Conti A., Zepf S.E., 1995, AJ 110, 1164 

Bohlin R.C., Savage B.D., Drake J.F., 1978, ApJ 224, 132 

Bruzual G.A., Chariot S., 1993, ApJ 405, 538 

Burstein D., Heiles C, 1984, ApJS 64, 601 

Carollo CM., Danziger I. J., 1994, MNRAS 270, 743 

Carollo CM., de Zeeuw P.T., van der Marel R.P., Danziger 

I. J., Qian E.E., 1995, ApJ 441, L25 
Conti P.S., Vacca W.D., 1994 ApJ 423, L97 
Davies R.L., Burstein D., Dressier A., Faber S.M., Lynden-Bell 

D., Terlevich R.J., Wegner C, 1987, ApJS 64, 581 
Djorgovsky S., Santiago B.X., 1992, ApJ 391, L85 
Fritze v. Alvesleben U., Burkert A., 1995, A&A 300, 58 
Geisler D., Piatti A. E., Claria J. J., Minniti D., 1995, AJ, 109, 

605 

Goudfrooij P., de Jong T., 1995, A&A 298, 784 

Goudfrooij P., Hansen L., Jorgensen H.E., N0rgaard-Nielson 

H.U., 1994a, A&AS 105, 341 
Goudfrooij P., Hansen L., Jorgensen H.E., N0rgaard-Nielsen 

H.U., de Jong T., van den Hoek L.B., 1994b, A&AS 104, 

179 

Gourgoulhon E., Chamaraux P., Fouque P., 1992, A&A 255, 
69 

Guhathakurta P., Knapp G.R., van Gorkom J.H., Kim D.-W., 
1989, in "The Interstellar Medium in External Galaxies", 



ed. D.J. Hollenbach & A.T. Harley, NASA Conf. Publica- 
tion 3084, p. 26 
Harris W.E., 1991, ARAA 29, 543 

Harris W.E., 1994, in "The Globular Cluster-Galaxy Con- 
nection" ASP Conf. Series Vol. 48, eds. G.H. Smith & 
J. P. Brodie 

Harris W.E., 1995, in "Stellar Populations", IAU Symp. 164, 
eds. P.C van dor Kruit & G. Gilmore, Dodrecht: Kluwer, 
p.85 

Harris W.E., 1996, electronically published catalog (McMaster 
University) 

Harris W.E., van den Bergh S., 1981, AJ 86, 1627 
Holtzman J. A., et al., 1992, AJ 103, 691 

Kim D., Jura M., Guhathakurta P., Knapp G.R., van Gorkom 

J.H., 1988, ApJ 330, 684 
Kissler-Patig M., Richtler T., Hilker M., 1996, A&A 308, 704 
Knapp G.R., Guhathakurta P., Kim D.-W., Jura M., 1989, 

ApJS 70, 239 

Kumai Y., Hashi Y., Fujimoto M., 1993, ApJ 416, 576 
Landolt A.U., 1992, AJ 104, 340 
Lima-Neto C, Combes F., 1995, A&A 294, 657 
Lutz D., 1991, A&A 245, 31 

McLaughlin D.E., Harris W.E., Hanes D.A., 1993, ApJ, 409, 
L45 

Minniti D., Alonso M.V., Goudfrooij P., Jablonka P., Meylan 

C, 1996, ApJ in press 
Poulain P., Nieto J.-L., 1994, A&AS 103, 573 
Richtler T., 1995, in "Reviews of Modern Astronomy", Vol. 8, 

eds. G. Klare, Springer, p. 163 
Schwcizer F., 1987, in "Structure and Dynamics of Elliptical 

Galaxies", IAU Symp. 127, ed. P.T. de Zeeuw, Dordrecht: 

Reidel, p. 109 
Schweizer F., Seitzer P., 1992, AJ 104, 1039 
Schweizer F., Seitzer P., 1993, ApJ 417, L29 
Schwcizer F., Seitzer P., Faber S.M., Burstein D., Dalle Ore 

CM., 1990, ApJ 364, L33 
Seeker J.A., 1992, AJ 104, 1472 

Smail I., Hogg D.W., Yan L., Cohen J.G., 1995, ApJ 449, L105 
Tully R.B., 1988, Nearby Galaxies Catalog, Cambridge Uni- 
versity Press 

Vaucouleurs de G., Vaucouleurs de A., Corwin H.G., Buta R.J., 
Paturel C, Fouque P., 1991, Third Ref. Catalogue of Bright 
Galaxies, Springer, New York 

White S.D.M., 1994, in: Lecture Series of the Les Houches Sum- 
mer scholl on 'Dark Matter and Cosmology' August 1993, 
ed. R.Schaeffer, North Holland, in press, also preprint MPA 
831 

Whitmore B.C., Schweizer F., Leitherer C, Borne K., Robert 

C, 1993, AJ 106, 1354 
Whitmore B.C., Schweizer F., 1995, AJ 109, 960 
Zepf S.E., Ashman K.M., 1993, MNRAS 264, 611 
Zepf S.E., Geisler D., Ashman K.M., 1993, ApJ 435, L117 
Zepf S.E., Ashman K.M., Geisler D., 1995, ApJ 443, 570 



This article was processed by the author using Springer- Verlag 
DT E X A&A style file L-AA version 3. 



